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Abstract 


The  maximin  algorithm  adapts  the  weights  of  an  antenna  array  to 
provide  simultaneous  interference  suppression  and  beamforming  in  a 
frequency-hopping  communication  system.  The  algorithm  is  derived  and 
a  digital  implementation  is  presented.  Simulation  experiments  illustrate 
the  ability  of  the  maximin  algorithm  to  rapidly  form  deep  nulls  in  the 
directions  of  interference  sources.  For  the  full  exploitation  of  an  adaptive 
array  for  wideband  frequency-hopping  communications,  some  form  of 
frequency  compensation  is  necessary.  The  strengths,  limitations,  and 
variations  of  three  methods  are  examined. 


Preface 


Recent  advances  in  digital  technology  and  the  increased  possibility  of 
Army  applications  of  frequency  hopping  have  rekindled  our  interest  in 
the  maximin  algorithm,  which  was  originally  proposed  and  studied  more 
than  a  decade  ago.  This  report  modifies  and  extends  the  original  pub¬ 
lished  research  results.  Along  with  many  minor  improvements,  there  are 
four  major  revisions: 

1.  The  theory  is  derived  in  terms  of  complex  envelopes  rather  than  analytic 
signals,  because  complex  envelopes  are  directly  extracted  in  many  mod¬ 
ern  digital  implementations. 

2.  The  use  of  a  one-sided  monitor  filter  rather  than  a  two-sided  band-reject 
filter  is  recommended. 

3.  A  new  convergence  analysis  of  this  highly  nonlinear  algorithm  is 
presented. 

4.  Far  more  extensive  simulation  experiments  indicate  the  potential  power 
and  the  limitations  of  the  maximin  algorithm. 
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1.  Introduction 


The  maximin  algorithm  is  an  adaptive-array  algorithm  that  suppresses 
interference  before  it  enters  the  demodulator  of  a  frequency-hopping 
communication  system  and  thereby  provides  a  spatial  processing  gain 
that  supplements  the  inherent  processing  gain  of  the  frequency-hopping 
system  [1-3].  The  algorithm  discriminates  between  the  desired  signal  and 
the  interference  on  the  basis  of  the  distinct  spectral  characteristics  of 
frequency-hopping  signals.  The  maximin  algorithm  is  so  named  because 
the  desired  signal  is  enhanced  and  the  interference  is  suppressed  simulta¬ 
neously.  The  maximin  algorithm  is  a  blind  adaptive  algorithm.  Blind  algo¬ 
rithms  do  not  require  training  sequences,  decision-directed  adaptation,  or 
knowledge  of  the  direction  of  the  desired  signal,  but  instead  rely  on  the 
basic  characteristics  of  the  desired-signal  waveform.  The  maximin  algo¬ 
rithm  requires  only  that  the  frequency-hopping  pattern  must  be  known 
by  the  receiver. 
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2.  Derivation 


The  basic  configuration  of  an  adaptive  array  for  frequency-hopping 
communications  is.  displayed  in  figure  1.  The  frequency-hopping  replica, 
which  is  produced  by  a  synchronized  frequency  synthesizer,  is  used  to 
dehop  the  arriving  desired  signals  [4],  thereby  removing  the  time  varia¬ 
tions  in  the  carrier  frequency.  The  sampled  values  of  the  complex  envel¬ 
opes  of  the  dehopped  signals  are  extracted  by  the  initial  processors  and 
applied  to  the  adaptive  filter,  which  executes  the  maximin  algorithm. 

The  maximin  algorithm  can  be  derived  by  the  method  of  steepest  descent 
[4]  with  the  signal-to-interference-plus-noise  ratio  (SINR)  as  the  perform¬ 
ance  criterion.  The  desired  signal  and  the  interference  are  assumed  to 
arrive  at  an  adaptive  array  of  N  antennas.  The  desired  signal,  interference 
signals,  and  thermal  noise  are  modeled  as  independent  zero-mean,  wide- 
sense-stationary  stochastic  processes.  Let  X(i)  denote  the  discrete-time 
vector  of  the  complex  envelopes  of  the  N  antenna  outputs  after  each  one 
has  been  dehopped,  filtered,  and  sampled.  The  index  i  denotes  the  sample 
number.  The  vector  X(z)  can  be  decomposed  as 

X(z)  =  s(t)  +  n(z)  ,  (1) 

where  s(z)  is  the  vector  of  desired-signal  complex  envelopes,  and  n(z)  is 
the  vector  of  interference  and  thermal-noise  complex  envelopes.  Let  W 
denote  a  vector  of  complex  weights.  Let  an  asterisk  denote  the  complex 
conjugate,  T  the  transpose,  and  H  the  conjugate  transpose.  The  complex 
system  output  is 

y(i)  =  X(f)  =  ys(z)  +  yn(i)  ,  (2) 

where 

ys{i)  =  WT  s(0,  yn(i)  =  WT  n(0  .  (3) 

The  power  in  the  complex  signal  is  defined  to  be  one-half  the  expected 
value  of  the  magnitude  squared.  The  factor  of  one-half  is  inserted  for 


Figure  1.  Basic 
configuration  of 
adaptive  array  for 
frequency-hopping 
communications. 


To  demodulator 
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consistency  with  the  convention  that  the  power  in  the  complex  envelope 
s(i)  is  equal  to  the  power  in  the  original  bandpass  signal  (e.g.,  see  Proakis 
[5]).  If  W  is  regarded  as  deterministic,  the  desired-signal  output  power  is 


-  1  whrssw  , 


(4) 


where  Rss  is  the  desired-signal  correlation  matrix, 

Rss  =  E[s*(i)  sT(i)]  ,  (5) 


and  E[x]  denotes  the  expected  value  of  x.  The  interference-plus-noise 
output  power  is 


p  =!whr  w 

1  n  2  nn  ' 


where  Rnn  is  the  interference-plus-noise  correlation  matrix 


R„„  =  E[n*(i)nT(0]  . 


(7) 


TheSINRis 


whr,w 

Pn  WHR  W 


(8) 


To  use  the  method  of  steepest  descent  in  deriving  an  adaptive  algorithm, 
it  is  necessary  to  calculate  the  gradient  of  the  SINR  with  respect  to  the 
weight  vector.  In  terms  of  its  real  part  and  its  imaginary  part  Wj,  a 
complex  weight  vector  is  defined  as 


W  =  W R-jVfj, 


(9) 


where  /  =  v'-l  .  Let  VWR  and  Vj^j  denote  the  gradients  with  respect  to  VfR 
and  Wjr,  respectively.  The  complex  gradient  with  respect  to  W  is  defined 
as 

Vw  =  ^WR~Fwi  ■  (10) 

A  direct  calculation  using  equations  (4)  and  (6)  gives 


Vw P  =  P 


Rssw 


(ii) 


Let  Ts  denote  the  time  interval  between  samples.  Let  m  denote  the  num¬ 
ber  of  samples  between  weight  iterations.  The  time  interval  between 
weight  iterations  is  m  Ts.  The  method  of  steepest  descent  for  discrete-time 
systems  gives  the  recursive  equation  for  the  weight  vector: 


W(*  +  1)  =  W (k)  +  |i0(Jt )  V(fc)  (12) 


where  k  denotes  the  weight  iteration  number,  fUo(k)  is  a  scalar  sequence 
that  controls  the  rate  of  change  in  the  weight  vector,  andV(fc)  is  an  esti¬ 
mate  of  the  gradient  in  equation  (11)  at  iteration  k.  The  ideal  steepest- 
descent  algorithm  results  if  V(fc)  is  set  equal  to  equation  (11)  with  W(fc)  in 
place  of  W.  However,  this  algorithm  requires  the  estimation  of  Rnn,  which 
is  unknown  in  general,  and  Rss,  which  is  unknown  in  the  absence  of 
information  about  the  direction  of  the  desired  signal.  To  avoid  the  direct 


3 


estimation  of  these  two  matrices,  several  approaches  are  possible.  If  the 
desired  signal  is  assumed  to  be  narrowband,  then  an  approximate  linear¬ 
ization  [4]  of  equation  (11)  leads  to  the  Howells- Applebaum  algorithm. 
Other  approximations  lead  to  the  Shor  feedback  loop  [6],  which  does  not 
completely  specify  an  adaptive  algorithm.  A  third  approach  leads  to  the 
maximin  algorithm,  which  is  developed  here. 

The  discrete-time  vector  s(z)  is  obtained  from  a  continuous-time  vector 
S(f)  of  complex  envelopes.  Each  component  of  S(f)  can  be  decomposed  as 

Si(t)  =  SiR(t)+jSi{t),  i  =  1, 2, ...  ,N  ,  (13) 

where  S,r(0  is  the  real  part  or  in-phase  component  of  S,-(f ),  and  S,-/(f)  is  the 
imaginary  part  or  quadrature  component  of  S,(f).  Each  S//(f),  i  =  1,2, , 
N,  is  a  delayed  version  of  the  first  component  Su(t),  which  is  modeled  as 
a  zero-mean,  wide-sense  stationary  process.  Two  fundamental  properties 
of  the  components  of  a  complex  envelope  are  [4,5] 

e[Si R(t)  S1R(t  +  r)]  =  E[sn(f)  su{t  +  t)]  (14) 

E[slR(t)  Sjj(f  +  r)j  =  -  EjSjj(t)  S1R(f  +  r)J  ,  (15) 

where  t  is  an  arbitrary  delay.  Using  these  equations,  we  obtain 

E[s(z)  sT(i)]  =  0  .  (16) 

The  desired-signal  component  of  the  adaptive-filter  output  can  be  decom¬ 
posed  as 

ys(0  =  ysR(0+;'ys/(0 ,  (17) 


where  ysn(i)  and  ys](i)  are  the  real  and  imaginary  parts  of  ys(i),  respec¬ 
tively.  If  the  weight  vector  is  treated  as  a  constant  W,  then  equations  (3), 
(16),  and  (17)  imply  that 


This  equation  and  equation  (5)  yield 

RssW  =  2E[s*(i)ysR(i)]  .  (19) 


Similarly,  if  the  interference  arriving  at  each  antenna  is  a  delayed  version 
of  the  interference  arriving  at  antenna  1,  and  if  the  thermal  noise  is  inde¬ 
pendent  in  each  array  branch,  then 

E[n(i)  nT(z)]  =  0  ,  (20) 

and,  hence. 


Run  W  =  2  E[n*(z)  ynR(j)]  , 


(21) 
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where  ynRp)  is  the  real  part  of  yn(i).  Substitution  of  equations  (19)  and  (21) 
into  equation  (11)  yields 


Vwp=2p 


E[s*3/sr]  £[n*y«K] 


where  the  time  index  has  been  omitted,  because  stationarity  implies  that 
the  expected  values  are  independent  of  time. 

To  derive  the  maximin  algorithm,  let  Ps(k)  and  Pn{k)  denote  estimates  of  Ps 
and  PM  at  iteration  k.  The  estimate  of  r  at  iteration  k  is 


<23) 

Let  s*ysR(k)  and  n*yn^(k)  denote  estimates  at  iteration  k  of  the  desired- 
signal  correlation  vector  E[s*i/sr]  and  the  interference-plus-noise 
correlation  vector  E[n *ynR],  respectively.  The  scalar  sequence  po(k)  in 
equation  (12)  ultimately  controls  the  stability  and  rate  of  convergence  of 
the  maximin  algorithm.  It  is  desirable  that  this  sequence  be  chosen  so  that 
as  p{k)  increases,  the  changes  in  W(fc)  will  be  small.  Thus,  we  set 

•  <24) 

where  a  is  a  constant.  Using  equations  (22)  and  (24)  and  the  preceding 
estimates  in  equation  (12),  we  obtain  the  basic  form  of  the  maximin 
algorithm: 

+  +  (25) 

where  k  >  0  and  W(0)  is  a  deterministic  initial  weight  vector.  The  subse¬ 
quent  convergence  analysis  and  simulation  results  confirm  that  the  choice 
given  by  equation  (24)  is  effective  and  robust,  provided  that  the  adaptation 
constant  a  is  appropriately  selected.  The  first  term  within  the  brackets  in 
equation  (25)  can  be  interpreted  as  a  signal  term  that  enables  the  algorithm 
to  direct  the  array  beam  toward  the  desired  signal.  The  second  term 
within  the  brackets  is  a  noise  term  that  enables  the  algorithm  to  null 
interference  signals. 

The  remaining  issue  is  the  choice  of  estimates  for  s*ysR(k) ,  n*y  nR(k) , 

Ps(k) ,  and  Pn(k) .  The  specific  nature  of  frequency-hopping  signals  allows 
these  estimates  to  be  made  without  depending  on  known  steering  vectors 
or  reference  signals,  as  shown  in  the  next  section. 
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3.  Implementation 

Figure  2  illustrates  the  principal  components  of  a  digital  implementation 
of  the  processing  behind  antenna  1  of  an  adaptive  array;  the  processing 
behind  the  other  antennas  is  similar.  The  front-end  devices  include  a 
bandpass  filter  that  excludes  noise  outside  the  hopping  band  and  a  low- 
noise  amplifier.  After  dehopping  from  a  variable  frequency  to  a  fixed 
intermediate  frequency  (IF),  the  received  signal  is  filtered  to  prevent 
aliasing  and  to  reject  noise.  The  signal  is  then  sampled  by  an  analog-to- 
digital  (A/D)  converter  at  a  rate  typically  four  times  the  intermediate 
frequency.  These  samples  are  applied  to  a  baseband  converter,  which 
typically  includes  a  Hilbert  transformer  and  a  complex  multiplier  [7],  to 
produce  sampled  values  of  the  complex  envelope  of  the  intermediate- 
frequency  signal.  The  discrete-time  output  of  the  converter  is  passed 
through  a  baseband  filter  with  a  bandwidth  equal  to  that  of  a  frequency 
channel.  The  filter  blocks  as  much  interference  and  noise  as  possible 
without  distorting  the  desired  signal.  Thus,  the  passband  of  the  frequency 
channel  being  used  by  the  frequency-hopping  desired  signal  is  translated 
to  baseband.  The  baseband-filter  output  is  the  first  complex  component  of 
X(i),  a  vector  applied  to  the  adaptive  filter. 

Figure  3  illustrates  the  digital  processing  required  by  the  maximin  algo¬ 
rithm.  If  the  interference  occupies  only  a  small  fraction  of  the  hopping 

x, 


"i 

Figure  2.  Dehopping  and  initial  processing. 

Figure  3.  Adaptive 
filter  that  executes 
maximin  algorithm. 
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band,  then  it  is  usually  absent  from  X(z)  and  x/r(i).  As  the  weights  con¬ 
verge  to  their  steady-state  values,  the  interference  component  of  i/r(0 
becomes  small,  even  if  the  interference  is  present  in  X(i).  Thus,  if  the  noise 
power  is  much  less  than  the  desired-signal  power,  then  a  suitable  esti¬ 
mate  of  the  desired-signal  correlation  vector  is 


*ys#)  “ 


km 

m  X  X*(z)  yR(z  , 

mi  =  (k-l)m  +  l  w  KW 


(26) 


where  is  the  real  part  of  the  adaptive-filter  output  y(z)  given  in 
equation  (2). 


An  analogous  estimate  for  n*ynR(fc)  is  not  possible,  because  when  the 
system  hops  into  a  part  of  the  band  with  interference  power,  then  both 
the  desired  signal  and  the  interference  are  imbedded  in  X(z).  Instead,  the 
receiver  observes  the  interference  and  noise  in  a  nearby  frequency  chan¬ 
nel  that  is  not  currently  being  used  by  the  desired  signal,  but  will  be  after 
a  subsequent  frequency  hop.  After  a  frequency  translation,  the  interfer¬ 
ence  and  noise  in  the  nearby  monitored  channel  are  extracted  by  each 
antenna's  monitor  filter  from  the  baseband  converter  output,  as  shown  in 
figure  2.  The  intermediate-frequency  filter  must  have  a  passband  large 
enough  to  encompass  both  the  downconverted  frequency  channel  being 
used  by  the  desired  signal  and  the  downconverted  monitored  channel. 


The  monitor  filter  has  a  passband  offset  by  f0  from  the  baseband  filter  and 
coinciding  with  a  downconverted  frequency  channel,  as  depicted  by  the 
transfer  function  sketched  in  figure  4.  The  parameter  /?  is  the  ratio  of  the 
bandwidth  of  the  monitor  filter  to  the  bandwidth  of  the  baseband  filter. 
Alternatively,  the  monitor  filter  may  be  a  band-reject  filter  with  two 
passbands  symmetrically  located  and  an  approximate  notch  over  the 
passband  of  the  baseband  filter  [1-3].  The  monitor  filter  allows  the 
adaptive  filter  to  monitor  the  interference  that  will  be  present  in  the 
baseband  filter  when  the  carrier  frequency  of  the  desired  signal  eventu¬ 
ally  coincides  with  that  of  the  interference.  Whenever  this  interference  is 
observed  at  the  output  of  the  monitor  filter,  the  adaptive  weights  are 
changed  to  cancel  the  interference.  There  is  negligible  spillover  of  the 
desired  signal  into  the  output  of  the  monitor  filter  if  the  filter  offset  is 
fo  =  2B  or  f0  =  -2B,  where  B  is  the  bandwidth  of  both  the  baseband  filter 
and  a  frequency  channel.  Subsequently,  the  single-sided  monitor  filter  of 


Figure  4.  Transfer 
functions  of  filters  for 
complex  envelopes. 


Magnitude 

Baseband  Monitor 
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figure  4  with  /?  =  1  is  assumed,  because  then  the  monitor  filter  passes 
almost  exactly  the  same  interference  waveforms  as  are  subsequently 
passed  by  the  baseband  filter. 


When  the  carrier  frequency  is  at  or  near  one  of  the  ends  of  the  hopping 
band,  the  monitor  filter  will  not  be  able  to  observe  interference  in  the 
hopping  band.  However,  this  effect  will  be  of  negligible  importance  when 
there  are  a  large  number  of  carrier  frequencies. 

Let  n(z')  denote  the  vector  of  N  discrete-time  outputs  of  the  monitor  filters 
and  define 


?»(>') = wTM  *(0  •  <27> 

The  vector 


nc(z)  =  n(f)exp(-;2it/0Ts)  (28) 

provides  an  estimate  of  the  baseband-filter  input  n(z)  that  will  occur  when 
the  desired-signal  carrier  frequency  is  such  that  the  interference  and  noise 
currently  in  the  monitor  filter  enter  the  baseband  filter.  Treating  the 
weight  vector  as  a  constant  so  that  equation  (21)  is  valid  and  assuming 
that  Rn„  ~  L[ne(z)  nj(z)j ,  we  obtain 


E[n10  3/ni?(0]  -  \  £[ne(z)  nJ(0]  w<0  • 

From  equation  (20),  it  is  reasonable  to  assume  that  E|ne(j)  n J(/)J  ~  0 
assumption  and  equations  (27)  to  (29)  yield  J 


(29) 

This 


E[n*(0y„R(0]-E 


n'O'K  \  nT(0  W(z)  +  i  hn(i)  W\i) 


1  aH/.-v 


(30) 


ft*(0  ?«r(0]  / 


where  y„R(i)  is  the  real  part  of  yn[i) .  This  equation  implies  that  a  suitable 
estimate  is 


-  km 


(31) 


The  estimate  Ps(k)  can  be  simplified  by  observing  from  a  straightforward 
calculation  using  equations  (3),  (16),  and  (17)  that  e[]/sk]  =  E[Vsl]  / so  that 


equation  (4)  implies  that  Ps  =  E[y^Rj «  E[y|] .  Similarly,  Pn  =  E[y^Rj 

" 2  1  •  Therefore,  suitable  power  estimates  are  obtained  by  computing 


*  E 


(32) 

km 

(33) 

In  applications  where  the  desired  signal  may  sometimes  be  absent,  it  may 
be  necessary  to  lower-bound  Ps(fc)  by  a  small  positive  number  chosen  to 
prevent  a  possible  singularity  in  the  maximin  algorithm. 

The  estimates  in  equations  (26),  (31),  (32),  and  (33)  complete  the  specifica¬ 
tion  of  the  maximin  algorithm  given  by  equation  (25).  If  the  underlying 
assumptions  are  changed,  then  different  versions  of  these  estimates,  and 
hence  of  the  maximin  algorithm,  result.  For  example,  if  the  interference 
occupies  large  spectral  regions  of  the  hopping  band,  then  usually  x/R(i) 
will  contain  significant  interference  power  before  the  weights  converge. 

To  compensate  for  this  power,  closely  related  alternative  versions  of  the 
maximin  algorithm  can  be  derived  [1-3].  However,  many  simulation 
experiments  indicate  that  the  version  of  the  algorithm  derived  above  is 
the  most  flexible,  nearly  always  gives  a  similar  or  better  performance  than 
alternative  versions,  and  is  computationally  the  simplest. 

During  acquisition  of  the  frequency-hopping  pattern,  the  frequency 
synthesizer  produces  different  frequencies  from  those  of  the  received 
frequency-hopping  pulses.  Consequently,  the  desired  signal  is  usually 
absent  from  the  outputs  of  the  baseband  filters,  and  the  maximin  algo¬ 
rithm  must  be  modified.  During  acquisition,  a  suitable  adaptive  algo¬ 
rithm  is  the  recursive  suppression  algorithm  [8],  which  tends  to  suppress 
signals  arriving  from  any  direction.  Cancellation  of  the  desired  signal  is 
avoided  because  of  its  infrequent  appearance  in  the  outputs  of  the 
baseband  and  monitor  filters. 

During  acquisition,  the  adaptive  system  forms  nulls  in  the  directions  of 
the  interference  signals,  thereby  assisting  the  receiver  in  synchronizing 
with  the  received  frequency-hopping  signal.  Although  a  grating  null 
might  be  inadvertently  placed  on  the  desired  signal  during  acquisition, 
the  probability  of  this  event  can  usually  be  kept  small  by  appropriate 
placement  of  the  antenna  elements  [6].  After  acquisition  is  confirmed  by 
the  receiver,  it  activates  the  maximin  algorithm. 
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4.  Convergence  Analysis 


A  frequency-hopping  signal  consisting  of  pulses  transmitted  at  many 
different  frequencies  is  a  wideband  signal.  However,  the  frequency 
channel  associated  with  a  single  pulse  is  usually  narrow.  Thus,  although  a 
frequency-hopping  signal  may  hop  over  a  wide  band,  it  has  a  narrow 
instantaneous  bandwidth.  For  a  narrowband  desired  signal,  the  vector  of 
sampled  complex  envelopes  is  well  approximated  by  [4,6] 

s(i)  =  8j(0  S0  ,  (34) 

where  Si(i)  is  the  sampled  complex  envelope  at  antenna  1,  which  serves  as 
a  reference  point,  and  So  is  a  steering  vector.  If  the  antenna  patterns  are 
identical  and  the  antennas  are  close  enough  that  the  signal  amplitudes  at 
all  the  antennas  are  nearly  the  same,  then  component  i  of  the  steering 
vector  is  exp(-/  2  nfc  Tj ),  where  fc is  the  center  frequency  and  Tj,  2  <i<  N, 
is  the  arrival-time  delay  at  antenna  i  relative  to  antenna  1.  After  the 
substitution  of  equations  (5)  and  (34)  into  equation  (8),  the  maximization 
of  the  SINR  yields  the  optimal  weight  vector  [4,6] 

W0=>IR™Si.  t35’ 

where  r\  is  an  arbitrary  constant,  and  the  superscript  -1  denotes  the 
inverse.  For  those  pulses  that  hop  over  a  sufficiently  small  band  within 
the  hopping  band,  the  optimal  weight  vector  is  well  approximated  by 
equation  (35),  with  Sq  evaluated  at  the  center  frequency  of  the  small  band. 


The  highly  nonlinear  nature  of  the  maximin  algorithm  precludes  a  com¬ 
pletely  rigorous  convergence  analysis.  However,  with  enough  approxima¬ 
tions,  the  convergence  of  the  mean  weight  vector  to  Wq  can  be  demon¬ 
strated  and  bounds  on  the  adaptation  constant  can  be  derived.  Assume 
that  enough  algorithm  iterations  have  occurred  that  the  power  estimates 
are  well  approximated  by  the  values  they  would  have  if  the  weight  vector 
were  equal  to  Wq.  The  desired-signal  power  at  an  antenna  is 


P  .  =  1  F 
1  si  2  C 


(36) 


Using  equations  (4),  (5),  (34),  (35),  and  (36),  and  the  fact  that  R„„  is  Hermi- 
tian,  we  obtain 


p#>  -  ■ V =  I  <  Rss  w0  =  |  n  I2  P,i  (si  R~l  s‘0f  .  (37) 

Similarly, 

-  Pm  =  \  w"  R„„  w0  =  M-  s\  r;’sJ>  .  (38) 

Equations  (37)  and  (38)  imply  that 


10 


Assuming  the  zero-mean  stationarity  of  the  signals  and  the  independence 
of  the  desired-signal  and  the  interference-plus-noise,  we  obtain  from 
equations  (26),  (1),  and  (2)  that 

£[i=^W]  =  E[x*(i)yK(i)]  (40) 

=  E[s»(i)ysj(i)] +  £[»*(')  ynn(')]  • 

Equation  (25)  indicates  that  W (k  +  1)  depends  on  S (l)  and  n(Z)  for  l<k,  but 
not  directly  on  S (k  +  1)  and  n (k  +  1).  Therefore,  we  make  the  approxima¬ 
tion  that  W (k)  is  statistically  independent  of  S (k)  and  n(k).  A  calculation 
similar  to  equation  (18)  using  equations  (16)  and  (20)  then  yields 

£[^«]-5(R-  +  R»»)£[W®]'  <41) 

Similarly,  equations  (30)  and  (31)  imply  that 

Taking  the  expected  value  of  both  sides  of  equation  (25),  substituting 
equations  (37)  to  (42),  and  simplifying  algebraically,  we  obtain 

E[W(fc  +  1)]-/i-  2^jr  [(p. - 1)  - *«]  j  E|W«]  ,  <«> 

where  k> 0  and  E[W(0)]  =  W(0).  Assuming  that  p0  »  1,  we  set p0-l~  p0 
and  define 

D  =  Po  Rn«  “  Rss  •  (44) 

We  obtain  the  approximate  result 

E[W(*  +  1)]  =  |l  -  D  J  E[W(*)]  .  (45) 

A  straightforward  calculation  using  equations  (5),  (34),  (36),  and  (39) 
yields 

D  Rnn  Sq  =  0  .  (46) 

“1  * 

This  equation  indicates  that  Rnn  S0  is  an  eigenvector  of  D,  and  the  corre¬ 
sponding  eigenvalue  is  0.  Since  D  is  Hermitian,  it  has  a  complete  set  of  N 
orthogonal  eigenvectors,  one  of  which  is  R^  SQ .  Thus,  we  can  make  the 
decomposition 

E[w(*)]  =  m  r;1  So  +  afk)  ef  ,  (47) 

where  each  a,(k)  and  t](k)  are  scalar  functions  and  e,-  is  one  of  the  N  - 1 
eigenvectors  orthogonal  to  R^  SQ  .  Substituting  equation  (47)  into  equa¬ 
tion  (45)  and  using  the  orthogonality  of  the  eigenvectors,  we  obtain 

r)(k  +  l)  =  rj(k)  =  T}  (48) 

«,<*  +  1)  =  f  1  -  Yfp~  Ai]  «,<*)/  i  =  1'2 . N ~ 1  '  (49) 
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where  rj  is  an  arbitrary  constant  and  A,-  is  the  eigenvalue  corresponding  to 
e*.  Equations  (48),  (47),  and  (35)  give 

E[m)]  =  WQ  +  Nz\l(k)ei  ,  (50) 

which  implies  that  E[W(/c)]  converges  to  Wo  if  and  only  if  each  a,(fc) 
converges  to  0.  From  equation  (49),  it  follows  that 


fl^)  =  (1-2^^)  k~°  • 

Therefore,  a,{fc)  converges  to  0  if  and  only  if 


1- 


a 


2  P0Pt 


SO 


<1  ,  i  =  l,2,...,N-l 


(51) 


(52) 


As  shown  subsequently,  each  A,-  is  positive.  Equation  (52)  yields  the 
necessary  and  sufficient  convergence  condition  for  the  convergence  of 
£[W(*)]: 


0<  a< 


*p0p: 


so 


max 


where  is  the  largest  eigenvalue  of  D. 


(53) 


It  remains  to  show  that  each  A,-  is  positive.  Since  the  SINR  is  maximized 
when  W  =  Wq,  it  follows  from  equation  (8)  that 


wHR„w 

W«R„„W 


5  P, 


(54) 


where  the  inequality  is  strict  when  W  *  W0.  Combining  equations  (44) 
and  (54)  gives 

WhDW>0,  (55) 


which  proves  that  the  Hermitian  matrix  D  is  nonnegative  definite.  Hence 
all  its  eigenvalues  must  be  nonnegative.  Since  only  W  =  Wo  gives  an 
equality  in  equation  (55),  setting  W  =  e,  implies  that 

ef* De,->0  ,  i  =  l,2,.  ..,N-1  .  (56) 

Therefore,  A,-  =  0  is  impossible  and 

A/>0  ,  i  =  1, 2, . . . ,  N  - 1  .  (57) 


The  sum  of  the  eigenvalues  of  a  square  martrix  is  equal  to  its  trace.  Thus, 
^max  -  Xj  =  tr(E>)  =  P0h(Rnn)  ~  ^(RSs) 

S Potr{Rnn )  ' 


(58) 


where  tr( A)  denotes  the  trace  of  A.  Consequently,  equation  (53)  implies 
that 


0  <a< 


4  P 


so 


fr  R 


nn 


(59) 
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is  a  sufficient  condition  (but  not  a  necessary  condition)  for  the  conver¬ 
gence  of  the  mean  weight  vector  to  the  optimal  weight  vector. 

Assume  the  interference-plus-noise  power  in  each  component  of  n(i)  is 
equal  to  Pni.  Then  tr(Rnn)  =  2 N  Pni.  Let  |  z  |  denote  the  Euclidean  norm  of  a 
vector  z.  From  equations  (5),  (34),  (36),  and  (37)  and  the  Schwarz  inequal¬ 
ity,  it  follows  that  Pso  <  |  W0 12  |  S0 1  Psi .  If  each  component  of  S0  is  a 
complex  exponential,  |  S0 1  =N  and  Pso  <  N  |  W0 1  Psi .  Substituting 
E[W(0)]  =  W(0)  into  equation  (50)  and  using  the  orthogonality  of  the 
eigenvectors  gives  |  W0 1  <  |  W(0)  | .  Therefore,  the  upper  bound  au  =  4  Pso/ 
tr  Rnn  in  equation  (59)  is  approximately  bounded  by 

«„<2|W(0)|2Pl  ,  (60) 

where  Pl  =  Psj/Pni  is  the  SINR  in  a  single  baseband-filter  output  prior  to 
the  adaptive  filtering  and  equals  the  SINR  per  frequency  channel  at  an 
antenna  input.  Since  the  bound  in  this  inequality  is  easily  calculated,  it 
provides  a  useful,  albeit  imperfect,  guide  in  selecting  the  appropriate 
value  of  the  adaptation  constant. 
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5.  Description  of  Simulation 

In  the  simulation,  the  antenna  array  consists  of  four  omnidirectional 
antennas  located  at  the  vertices  of  a  square.  The  symmetry  of  this  archi¬ 
tecture  ensures  a  relatively  mild  dependence  of  the  array  response  on  the 
signal  directions  and  allows  full  azimuthal  coverage.  The  edge  length  is 
equal  to  1  or  2  times  the  wavelength  X  corresponding  to  the  center  fre¬ 
quency  of  the  desired  signal,  3  GHz.  Each  interference  source  is  in  the 
plane  of  the  array.  The  desired  signal  arrives  from  a  direction  perpendicu¬ 
lar  to  one  of  the  edges,  and  this  direction  is  defined  to  be  0°.  All  signals 
are  assumed  to  arrive  as  plane  waves.  The  system  parameters  are  summa¬ 
rized  in  table  1. 

The  frequency-hopping  signal  has  a  randomly  chosen  carrier  frequency 
within  the  hopping  band  and  is  modulated  by  binary  minimum-shift 
keying  (MSK).  The  sequence  of  data  bits  is  randomly  generated  at  the  rate 
of  100  kbps.  The  hop  dwell  time  is  1  ms.  The  thermal  noise  at  the  output 
of  each  intermediate-frequency  filter  following  an  antenna  element  is 
modeled  as  filtered  white  Gaussian  noise.  The  ratio  of  the  signal  power  to 
the  thermal-noise  power  is  defined  relative  to  a  single  frequency  channel 
and  is  set  equal  to  20  dB.  Since  the  frequency-hopping  signal  is  assumed 
to  be  noncoherent  in  carrier  phase  from  one  hop  to  the  next,  each 
frequency-hopping  pulse  has  an  initial  phase  that  is  uniformly  distrib¬ 
uted  over  the  interval  \-n,  n\.  The  hopping  frequencies  are  separated  by 
100  kHz  and  spread  uniformly  over  the  total  hopping  band,  which  occu¬ 
pies  either  30  or  300  MHz.  Thus,  there  are  either  300  or  3,000  contiguous 
frequency  channels.  The  total  interference  power  due  to  all  interference 
signals  is  equal  to  10  times  the  number  of  frequency  channels,  which 
maintains  a  constant  value  of  total  interference  power  per  frequency 
channel.  Thus,  pi  is  slightly  less  than  -10  dB.  Each  interference  signal  has 
the  same  power  and  distributes  tones  in  some  of  the  frequency  channels. 


Table  1.  System 
parameters. 


Parameter 

Array  antennas 
Array  edge  length 
Center  frequency 
Hop  dwell  time 
Data  rate 

Frequency  modulation 
Signal-to-noise  ratio 
Hopping  bandwidth 
Number  of  frequency  channels 
Monitor  filter  offset 
Sampling  rate 
Weight  iterations  per  hop 
Total  interference-to-signal  ratio 
Interference  type 
Number  of  hops  per  experiment 


Value 

4,  omnidirectional,  at  vertices  of  square 

1  or  2  wavelengths 

3  GHz 

1  ms 

100  kbps 

MSK 

20  dB  per  antenna  and  channel 

30  or  300  MHz 

M  =  300  or  3000 

200  kHz 

800  ksamples/s 

8 

10  M 

Tones  in  channels 
50 
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Perfect  synchronization  between  the  frequency-hopping  signals  at  all  the 
antenna  outputs  and  the  frequency  synthesizer  in  the  receiver  is  assumed. 
The  baseband  and  monitor  filters  are  modeled  as  digital  8-pole  Butter- 
worth  filters  with  the  3-dB  bandwidths  equal  to  100  kHz,  the  bandwidth 
of  a  frequency  channel.  The  monitor  filter  has  a  single  passband,  as 
depicted  in  figure  4,  with  parameters  f0  =  200  kHz  and  /3  =  1.  The 
sampling  rate  of  the  analog-to-digital  converter  is  800  kilosamples  per 
second.  There  are  8  weight  iterations  per  hop. 

In  each  of  the  subsequent  simulation  experiments,  the  initial  weight 
vector  is  W(0)  =  Wbeam/  where 

Wbeam  =  n  111]-  («) 

This  weight  vector  points  the  array  beam  or  gain  pattern  in  the  direction 
of  the  desired  signal.  Thus,  the  subsequent  simulation  experiments 
indicate  the  performance  of  the  maximin  algorithm  after  the  sudden 
appearance  of  interference  from  various  directions.  An  alternative  is  to  set 
W(0)  =  Womni,  where 

(62) 

This  weight  vector  corresponds  to  an  omnidirectional  gain  pattern.  Its  use 
as  the  initial  weight  vector  is  appropriate  if  one  wishes  to  assess  the 
convergence  of  the  maximin  algorithm  from  a  cold  start.  Since  the  algo¬ 
rithm  then  has  to  form  a  beam  in  the  direction  of  the  desired  signal  while 
simultaneously  attempting  to  null  the  interference  signals,  the  conver¬ 
gence  is  slowed  relative  to  what  it  would  be  if  W(0)  =  Wbeam-  Simulation 
experiments  verify  that  the  steady-state  performance  is  not  affected  by 
the  choice  of  Wbeam  or  Womni  as  the  initial  weight  vector,  so  Wbeam  is 
always  used. 

Inequalities  (59)  and  (60)  indicate  that  0  <a<au,  where  au  <  0.8  is  suffi¬ 
cient  for  convergence  if  the  interference  signals  approximate  stationary 
processes.  Simulation  experiments  confirm  that  the  maximin  algorithm 
does  converge  to  a  steady  state  when  0  <  a  <  0.8.  It  is  found  that  a  -  0.2 
usually  provides  close  to  the  fastest  convergence  and  the  best  perform¬ 
ance  when  each  interference  signal  comprises  tones  in  every  frequency 
channel. 

An  example  of  the  results  of  a  typical  simulation  trial  is  illustrated  in 
figure  5(a).  The  SINR  at  each  weight  iteration  is  calculated  by  averaging 
the  SINR  over  all  sampling  intervals  from  after  the  preceding  weight 
iteration  to  just  before  the  current  weight  iteration.  The  system  param¬ 
eters  are  array  edge  =  A  and  hopping  bandwidth  =  30  MHz.  Two  equal- 
power  interference  signals  arrive  from  40°  and  70°  relative  to  the  desired- 
signal  direction.  Each  of  these  signals  distributes  tones  in  all  the 
frequency  channels,  and  the  adaptation  constant  is  a  =  0.2.  The  SINR 
gradually  increases  until  the  adaptive  system  reaches  a  nearly  steady- 
state  condition  during  which  the  SINR  fluctuates  by  a  few  decibels. 
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Figure  5.  (a)  SINR  for 
single  simulation 
trial  with  two 
interference  signals, 
each  comprising 
tones  in  every 
frequency  channel, 
and  hopping 
bandwidth 
=  30  MHz,  (b)  array 
gain  pattern  at  end  of 
trial,  and  (c)  average 
SINR  for  20  trials. 


(a) 


Figure  5(b)  depicts  the  array  gain  pattern  at  the  end  of  the  simulation 
trial.  There  are  nulls  of  less  than  -20  dB  in  the  directions  of  the  interfer¬ 
ence  sources.  Figure  5(c)  shows  the  average  SINR,  which  is  defined  as  the 
SINR  obtained  by  averaging  over  20  simulation  trials  at  each  weight 
iteration. 

If  each  interference  signal  independently  distributes  tones  of  equal  power 
in  every  fourth  frequency  channel  to  cover  25  percent  of  the  hopping 
band,  then  the  simulation  results  are  quite  different.  As  shown  in  figures 
6(a)  and  (c),  the  convergence  is  slowed  and  the  fluctuations  in  SINR  are 
increased,  because  the  monitor  filter  only  observes  each  interference 
signal  25  percent  of  the  time.  Between  these  observations,  the  weights 
tend  to  slowly  drift  toward  the  values  they  would  have  without  the 
interference.  This  drift  and  the  occasional  large  interference  levels  force 
the  lowering  of  the  adaptation  constant  to  obtain  good  algorithm 
convergence,  and  a  =  0.05  is  used  to  generate  figure  6.  Figure  6(b) 
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Figure  6.  (a)  SINK  for 
single  simulation 
trial  with  two 
interference  signals, 
each  comprising 
tones  in  25  percent  of 
frequency  channels, 
and  hopping 
bandwidth 
=  30  MHz,  (b)  array 
gain  pattern  at  end  of 
trial,  and  (c)  average 
SINK  for  20  trials. 


Number  of  weight  updates 


indicates  that  the  array  gain  pattern  is  very  similar  to  that  of  the  preced¬ 
ing  example  as  a  nearly  steady-state  condition  is  approached. 


In  general,  interference  that  occupies  only  a  small  part  of  the  hopping 
band,  or  even  frequency-hopping  interference  signals,  can  be  suppressed 
by  the  maximim  algorithm  supplemented  by  an  error-correcting  code.  If 
the  interference  passes  through  the  monitor  filter  often  enough,  it  is 
suppressed  by  the  maximim  algorithm;  if  it  passes  through  only  occasion¬ 
ally,  then  it  is  accommodated  by  the  error-correcting  code. 


The  results  of  seven  representative  simulation  experiments  with  the  edge 
length  =  A  are  summarized  in  table  2.  The  interference  signals  distribute 
tones  in  all  the  frequency  channels  and  a  =  0.2  is  used.  Each  experiment 
comprises  20  trials  with  50  hops  and  400  weight  iterations  per  trial.  The 
first  column  gives  the  hopping  bandwidth.  The  second  column  gives  the 
arrival  angles  of  1, 2,  or  3  interference  signals  relative  to  the  desired-signal 
direction.  The  SINRs  for  the  last  40  weight  iterations  of  all  the  trials  are 
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averaged  to  obtain  the  final  SINR  and  are  used  to  calculate  the  standard 
deviation  of  the  final  SINR,  which  are  listed  in  the  third  and  fourth 
columns  of  the  table.  The  final  column  gives  the  crossing  number,  which  is 
defined  as  the  number  of  weight  iterations  required  for  the  average  SINR 
over  all  the  trials  to  exceed  a  threshold  value  equal  to  3  dB  less  than  the 
final  SINR.  The  crossing  number  provides  a  rough  measure  of  the  relative 
time  required  for  convergence  to  the  steady  state. 

The  first  row  of  table  2  shows  the  system  response  to  a  single  interference 
signal.  The  second  row,  which  is  obtained  from  the  same  experiment  that 
gave  figure  5,  indicates  that  there  is  slower  convergence  and  a  2.5-dB  loss 
in  final  SINR  due  to  the  presence  of  a  second  equal-power  interference 
signal.  The  third  and  fourth  rows  illustrate  the  further  losses  that  occur 
when  three  equal-power  interference  signals  are  present.  In  the  third  row, 
the  third  interference  signal  at  an  angle  -60°  fortuitously  arrives  near  a 
grating  null  naturally  generated  when  the  first  two  interference  signals 
are  nulled,  as  observed  in  figure  5(b).  In  the  fourth  row,  the  final  SINR  is 
much  smaller  because  the  third  interference  signal  arrives  in  the  natural 
grating  lobe.  The  fifth  and  sixth  rows  illustrate  the  decrease  in  the  final 
SINR  and  the  sharp  increase  in  the  standard  deviation  that  result  when 
the  bandwidth  increases  by  a  factor  of  10  while  the  total  interference 
power  remains  the  same.  The  final  row  illustrates  the  limitations  imposed 
by  the  resolution  of  the  array. 

The  resolution  of  an  adaptive  array  is  the  minimum  angular  separation 
between  the  interference  and  desired-signal  sources  that  can  be  accommo¬ 
dated  without  significant  performance  degradation.  The  resolution  is 
limited  by  the  angular  extent  of  a  null  formed  against  an  interference 
signal,  and  is  proportional  to  A/D,  where  D  is  the  aperture  or  separation 
of  the  antennas.  Thus,  to  decrease  the  resolution,  the  edge  length  of  the 
square  array  can  be  increased.  However,  the  increase  in  the  separation 
between  antennas  inevitably  causes  an  increase  in  the  number  of  grating 
nulls.  The  tendency  of  the  array  to  form  grating  nulls  at  certain  angles 
may  impede  the  formation  of  nulls  against  two  or  more  interference 
signals.  Furthermore,  since  a  decrease  in  the  resolution  causes  a  decrease 
in  the  angular  width  of  a  null,  there  is  a  decrease  in  the  fractional  band¬ 
width  that  can  be  accommodated  by  the  array,  as  explained  in  the  next 
section. 


liable  2.  Simulation 
results  with  edge 
length  =  A. 

Bandwidth 

(MHz) 

Direction  of 
interference  (°) 

Final 

SINR  (dB) 

Standard 
dev.  (dB) 

Crossing 

number 

30 

40 

25.7 

1.5 

5 

30 

40,70 

23.2 

1.3 

24 

30 

40, 70,  -60 

19.8 

1.9 

88 

30 

40, 70,  -50 

15.5 

2.1 

207 

300 

40 

23.7 

5.2 

6 

300 

40,70 

20.4 

4.3 

22 

30 

10 

19.0 

1.3 

128 

18 


Table  3  indicates  the  effects  of  an  increase  in  the  array  aperture  or  edge 
length  of  the  square  array.  The  superior  resolution  of  the  larger  array 
yields  greatly  improved  performance  when  the  interference  is  from 
approximately  the  same  direction  as  the  desired  signal.  However,  as 
indicated  by  the  sixth  row  in  the  table,  the  performance  advantage  usu¬ 
ally  disappears  when  the  hopping  bandwidth  is  large  and  there  are  two 
or  more  equal-power  interference  signals. 


Table  3.  Simulation 
results  with  edge 
length  =  2A. 


Bandwidth 

(MHz) 

Direction  of 
interference  (°) 

Final 

SINR  (dB) 

Standard 
dev.  (dB) 

Crossing 

number 

30 

40 

25.3 

1.5 

2 

30 

40,70 

24.5 

2.6 

15 

30 

40,  70,  -60 

19.9 

2.0 

58 

30 

40,  70,  -50 

23.6 

2.7 

22 

300 

40 

23.8 

2.2 

3 

300 

40,70 

14.0 

8.6 

12 

30 

10 

23.9 

1.3 

9 
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6.  Frequency  Compensation 


Adaptive  processing  may  be  significantly  impaired  if  the  fractional  band¬ 
width,  which  is  defined  as  the  hopping  bandwidth  divided  by  the  center 
frequency  of  the  hopping  band,  exceeds  a  few  percent  and  no  frequency 
compensation  is  used.  As  an  example,  the  hopping  bandwidth  in  the 
previous  example  of  figure  5  is  expanded  to  300  MHz,  a  10  percent  frac¬ 
tional  bandwidth,  but  all  the  other  conditions  remain  the  same.  The 
simulation  results  are  shown  in  figure  7.  The  significantly  increased  SINR 
fluctuations  are  caused  by  the  large  changes  in  the  carrier  frequency  of 
the  frequency-hopping  desired  signal  that  sometimes  occur.  Not  only  is 
convergence  slowed,  but  also  the  steady-state  SINR  is  reduced  relative  to 
the  results  in  figure  5. 


Frequency  compensation  is  sometimes  needed  because  the  antenna 
elements  produce  only  relative  phase  information.  Consider  two  antenna 


Figure  7.  (a)  SINR  for 
single  simulation 
trial  with  two 
interference  signals, 
each  comprising 
tones  in  every 
frequency  channel, 
and  hopping 
bandwidth 
=  300  MHz,  (b)  array 
gain  pattern  at  end  of 
trial,  and  (c)  average 
SINR  for  20  trials. 


(a) 
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Figure  8.  Two 
antenna  elements  of 
adaptive  array 
receiving  plane 
wave. 


elements  of  an  adaptive  array  receiving  a  plane  wave,  as  illustrated  in 
figure  8.  The  phase  shift  of  a  narrowband  signal  at  one  element  relative  to 
the  phase  of  the  signal  at  the  other  is 

6  =  f  D  sin  <f>  ,  (63) 

where  w  is  the  carrier  frequency  in  radians  per  second,  c  is  the  velocity  of 
an  electromagnetic  wave,  <f)  is  the  arrival  angle  of  the  plane  wave  relative 
to  a  line  perpendicular  to  the  straight  line  through  the  two  elements,  and 
D  is  the  distance  between  the  elements.  If  the  carrier  frequency  of  the 
received  signal  changes  from  ft)  to  ft)  +  A  ft),  then  the  phase  shift  changes  by 

(A0)1  =  ^  D  sin  0  .  (64) 

If  the  frequency  is  constant,  but  the  arrival  angle  changes  from  <j>  to 
<j)  +  A0,  then  the  phase-shift  change  is 

(A0)2  =  ®  D  sin  +  A^)  -  ®  D  sin  <j)  .  (615) 

If  the  adaptation  is  realized  by  the  use  of  a  single  adaptive  weight  behind 
each  antenna  element,  the  adaptive  system  cannot  distinguish  between 
the  two  possible  causes  for  a  change  in  the  phase  shift.  Thus,  a  change  in 
frequency  due  to  frequency  hopping  is  processed  as  if  a  change  occurred 
in  the  arrival  angle. 

Setting  (A0)!  equal  to  (A0)2,  we  conclude  that  the  equivalent  changes  in 
frequency  and  arrival  angle  are  related  by 

sin  +  A 0)  =  +  l]  sin  <j>  .  (66) 

This  equation  can  be  solved  for  A 0  as  a  function  of  A co  if  |  sin  </>  |  <  <0/ 

(ft)  +  A  co).  For  small  or  moderate  values  of  A  <p,  a  Taylor-series  expansion  of 
equation  (66)  yields 

A <p  =  tan  0,  |  A0  tan  <j)  |  <  <  1  .  (67) 

Suppose  that  an  adaptive  system  creates  an  approximate  spatial  null  in 
the  direction  of  a  source  of  interference  at  the  instantaneous  hopping 
frequency  ft) .  If  the  instantaneous  frequency  of  the  desired  signal  changes 
to  (o  +  Aft>  and  the  interference  source  produces  significant  power  at 
0  +  A  ft),  then  the  equivalent  A<j>  for  the  interference  is  given  by 

Plane  wave 
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equations  (66)  and  (67).  If  A  co  is  sufficiently  large,  A <j)  may  be  larger  than 
the  angular  width  of  the  original  null.  Consequently,  after  the  carrier 
frequency  of  the  desired  signal  hops  to  co  +  A®,  the  interference  is  not 
immediately  nulled,  but  further  adaptation  is  required  to  again  establish 
a  spatial  null.  Thus,  if  the  fractional  bandwidth  is  large,  then  unless 
(1)  the  angular  width  of  a  null  is  large,  (2)  A®  is  restricted  to  small  values, 
or  (3)  some  type  of  frequency  compensation  is  adopted,  the  benefit  of  the 
adaptive  processing  may  diminish  significantly. 

The  angular  width  of  a  null  can  be  increased  by  decreasing  the  separation 
among  antennas  in  the  array,  but  then  the  resolution  may  be  increased  to 
an  unsatisfactory  level.  The  restriction  of  the  frequency  change  after  a  hop 
A(o  to  a  small  value  keeps  the  interference  near  the  center  of  a  spatial  null 
and  makes  it  possible  for  the  adaptive  system  to  rapidly  increase  the 
interference  rejection.  However,  the  variety  of  hopping  patterns  is  re¬ 
duced,  which  diminishes  the  resistance  to  frequency-selective  fading  and 
increases  the  susceptibility  to  wideband  interference  or  jamming.  Thus, 
frequency  compensation  is  the  preferred  alternative. 

Three  basic  methods  of  frequency  compensation  for  frequency  hopping 
are  parameter-dependent,  spectral,  and  anticipative  processing. 
Parameter-dependent  processing  introduces  additional  parameters  in  an 
attempt  to  increase  the  angular  extent  and  spectral  bandwidth  of  the 
nulls.  For  the  maximin  algorithm,  a  compatible  realization  applies  v 
successive  discrete-time  outputs  of  each  baseband  filter  and  each  monitor 
filter  as  inputs  to  the  adaptive  filter,  thereby  increasing  the  number  of 
weights  in  figure  3  by  a  factor  of  v.  However,  the  filters  sharply  limit  the 
degree  of  spectral  or  beam  shaping  that  might  be  possible  using  the 
successive  outputs.  Simulation  results  confirm  that  parameter-dependent 
processing  is  ineffective  for  the  maximin  algorithm. 

Spectral  processing  is  based  on  dividing  the  hopping  band  into  a  number 
of  spectral  regions  and  adapting  independently  when  the  carrier  fre¬ 
quency  is  in  one  of  the  regions.  The  weights  associated  with  each  spectral 
region  are  stored  in  a  memory  within  the  weight  processor  of  figure  3.  At 
the  end  of  the  signal  dwell  time  at  a  specific  frequency,  the  weights  of  the 
adaptive  filter  are  transferred  to  the  memory.  Then  the  weights  associated 
with  the  new  carrier  frequency  are  transferred  from  the  memory  to  the 
filter.  These  weights  are  updated  during  the  time  interval  at  the  new 
frequency.  Transfers  between  the  memory  and  the  filter  are  controlled  by 
the  frequency-hopping  code  generator.  The  number  of  spectral  regions  Ns 
is  ultimately  limited  by  the  proportional  increase  in  the  required  number 
of  iterations  to  converge  to  steady  state,  because  only  1/NS  of  the 
frequency-hopping  pulses  are  associated  with  each  region.  However,  this 
increase  is  slowed  by  the  contraction  of  each  region  as  Ns  increases,  which 
expedites  the  convergence  within  each  region.  The  relatively  slow  conver¬ 
gence  rate  of  spectral  processing  is  its  major  disadvantage;  its  major 
advantage  is  its  relatively  simple  implementation. 


An  anticipative  adaptive  system  begins  adaptation  toward  the  optimal 
weights  for  a  carrier  frequency  before  that  frequency  is  transmitted. 
Anticipative  processing  with  an  auxiliary  processor  operating  in  parallel 
with  the  main  processor  is  illustrated  by  figure  9.  A  time-advanced 
frequency-hopping  replica  hops  one  hopping  period  ahead  of  the  replica 
used  for  dehopping  the  desired  signal.  While  the  main  adaptive  filter 
produces  the  output,  the  auxiliary  filter  adapts  the  weights  corresponding 
to  the  next  hopping  frequency.  After  each  hop,  the  weights  associated 
with  the  new  hopping  frequency  are  transferred  from  the  auxiliary  filter 
to  the  main  filter.  Transfers  are  triggered  by  the  clock  that  controls  the 
frequency-hopping  carrier  transitions.  The  weight  vector  in  the  auxiliary 
filter  is  not  necessarily  updated  at  the  same  times  as  the  weight  vector  in 
the  main  filter.  The  auxiliary  filter  may  update  its  weight  vector  during 
the  switching  times  between  hops,  but  in  the  simulation,  the  duration  of 
the  switching  time  is  assumed  to  be  negligible. 

The  desired  signal  is  not  present  in  the  auxiliary  processors  because  they 
monitor  the  frequency  channel  of  the  next  hop.  Consequently,  baseband 
filters  can  be  used  as  monitor  filters  to  extract  various  estimates  of  the 
interference-plus-noise.  These  baseband  filters  allow  narrower 
intermediate-frequency  filters  and,  hence,  a  lower  sampling  rate  in  the 
auxiliary  processor  than  in  the  main  processor. 


Let  W a(l)  denote  the  anticipative  weight  vector  computed  by  the  auxiliary 
adaptive  filter  at  iteration  l.  The  recursive  equation  for  Wfl(Z  +  1)  is  de¬ 
rived  by  modifying  the  maximin  algorithm  to  allow  for  the  absence  of  the 
desired  signal  in  the  auxiliary  filter  and  for  the  measurements  by  the 
auxiliary  filter.  Thus,  the  signal  term  within  the  brackets  in  equation  (25) 
is  discarded,  and  the  remaining  equation  is  simplified.  The  anticipative 
maximin  algorithm  is 


wa(/  +  i)  =  wfl(/)- 


m 


(68) 


where  a\  is  the  anticipative  adaptation  constant  and  the  subscript  a  on  the 
bracketed  factor  indicates  that  it  is  measured  by  the  auxiliary  adaptive 


Figure  9.  Basic 
configuration  of 
adaptive  array  with 
anticipative 
processing. 


To 

demodulator 
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filter,  which  uses  the  output  of  a  frequency  channel  that  will  not  be 
observed  by  the  main  adaptive  filter  until  the  next  hop.  The  desired- 
signal  power  estimate  PJl)  is  the  most  recent  estimate  computed  by  the 
main  adaptive  filter  and  is  assumed  to  be  nearly  independent  of  the 
frequency  channel.  The  bracketed  factor  is  defined  by  equation  (31),  with 
fi(0  interpreted  as  the  discrete-time  vector  of  complex  envelopes  at  the 
input  of  the  auxiliary  adaptive  filter  and  ynR{i)  as  the  real  part  of  the  filter 
output. 


The  weight  vector  in  the  main  adaptive  filter  is  updated  by  computing 
equation  (25),  except  at  sampling  instants  occurring  during  switching 
times.  At  these  instants,  the  weight  vector  in  the  main  adaptive  filter  is  set 
equal  to  that  of  the  auxiliary  adaptive  filter.  The  switching  times  occur 
when  k  =  nk0  in  the  main  filter  and  i  =  nk\  in  the  auxiliary  filter,  where  ko 
is  the  number  of  iterations  per  hop  in  the  main  filter,  k\  is  the  number  of 
iterations  per  hop  in  the  auxiliary  filter,  and  n  =  1, 2, ...  is  the  hop  num¬ 
ber.  Thus,  the  maximin  algorithm  in  the  main  filter  becomes 


W(fc+l)  =  W(jt) 


+  SL. 

P{k) 


SX#) 

w 


nX#) 

KM 


WK)  =  W«(nfc  l)'  n=0'1 - 


k  +  ljt  nkQ 


(69) 


The  anticipative  maximin  algorithm,  which  is  closely  related  to  the 
recursive  suppression  algorithm  [8]  for  direction  finding,  attempts  to  null 
all  signals  with  power  in  the  hopping  band,  but  does  not  null  the  desired 
signal  because  it  is  never  present  in  the  baseband  filters  of  the  auxiliary 
adaptive  filter.  This  absence  ultimately  is  due  to  the  system's  knowledge 
of  the  frequency-hopping  pattern. 


One  might  consider  a  version  of  the  anticipative  maximin  algorithm  in 
which  the  signal  term  in  equation  (25)  is  retained.  The  signal  term  cannot 
be  computed  by  the  auxiliary  adaptive  filter  because  it  does  not  receive 
the  desired  signal.  Therefore,  the  signal  term  computed  by  the  main 
adaptive  filter  is  used  to  enable  beamforming  by  the  auxiliary  adaptive 
filter  in  the  direction  of  the  desired  signal.  However,  the  signal  term  does 
not  include  frequency  compensation,  and  any  beamforming  tends  to 
erode  the  nulls  when  interference  is  absent  in  the  monitor  filter  outputs. 
Simulation  results  indicate  that  this  algorithm  version  can  accommodate 
larger  adaptation  constants  and,  hence,  expedites  the  convergence  when 
the  interference  signals  occupy  all  the  frequency  channels,  but  its 
performance  is  inferior  when  interference  signals  occupy  only  a  fraction 
of  the  hopping  band.  Thus,  equation  (68)  is  the  preferred  form  of  the 
anticipative  maximin  algorithm. 


The  results  of  nine  simulation  experiments  are  displayed  in  table  4,  where 
the  edge  length  is  A,  the  hopping  bandwidth  is  300  MHz,  and  the  frac¬ 
tional  bandwidth  is  10  percent.  The  first  three  rows  in  table  4  give  results 
for  no  frequency  compensation  and  a  =  0.2  as  a  reference  for  comparison. 
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Table  4.  Simulation 
results  with  edge 
length  =  A  and  total 
hopping  bandwidth 
=  300  MHz. 


The  next  three  rows  indicate  that  spectral  processing  with  three  regions 
and  a  =  0.2  provides  an  increased  final  SINR  at  the  cost  of  a  large  increase 
in  the  crossing  number.  The  standard  deviation  is  substantially  reduced 
for  one  or  two  interference  signals.  Using  two  regions  is  slightly  less 
effective.  Further  gain  in  the  final  SINR  and  a  reduction  in  the  standard 
deviation  occurs  as  the  number  of  hops  is  increased  beyond  the  50  used 
in  generating  table  4,  because  the  weights  of  each  of  the  three  spectral 
regions  approach  steady-state  values.  This  additional  gain  increases  with 
the  number  of  interference  signals. 

In  the  simulation  experiments  with  anticipative  processing,  the  initial 
value  of  the  auxiliary  weight  vector  is  Wfl  (0)  =  [111  1],  and  the  number 
of  iterations  per  hop  isk1  =  k0  =  8.  For  one  or  two  interference  signals,  the 
adaptation  constants  are  a  -  cc\  =  0.08;  for  three  interference  signals, 
cc—  cc-y  =  0.02.  Many  other  combinations  of  k\,  oc,  and  <X\  were  tried,  but 
none  gave  an  improved  overall  performance  compared  with  the  results 
listed  in  the  last  three  rows  of  table  4.  The  results  indicate  that  anticipative 
processing  provides  better  performance  in  both  the  final  SINR  and  the 
standard  deviation  than  spectral  processing  or  no  compensation, 
although  the  convergence  rate  of  anticipative  processing  against  three 
interference  signals  is  very  slow. 

The  simulation  results  indicate  that  when  the  fractional  bandwidth  is 
10  percent,  either  spectral  or  anticipative  processing  usually  provides  a 
significant  benefit.  Anticipative  processing  has  the  advantage  of  usually 
providing  the  best  performance,  but  spectral  processing  is  simpler  to 
implement. 


Compensation 

method 

Type 

Direction  of 
interference  (°) 

Final 

SINR  (dB) 

Standard 
dev.  (dB) 

Crossing 

number 

None 

— 

40 

23.7 

5.2 

6 

None 

— 

40,70 

20.4 

4.3 

22 

None 

— 

40, 70,  -60 

13.6 

5.6 

63 

Spectral 

Ns- 3 

40 

25.1 

2.7 

6 

Spectral 

Ns  =  3 

40,70 

22.1 

2.4 

70 

Spectral 

fcj  =  8 

40, 70,  -60 

14.9 

5.8 

142 

Anticipative 

fc1  =  8 

40 

26.0 

1.3 

12 

Anticipative 

k,  =  8 

40,70 

23.4 

1.4 

32 

Anticipative 

JV 

It 

00 

40, 70,  -60 

16.0 

5.1 

190 
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The  maximin  algorithm  is  a  blind  adaptive  algorithm  that  simultaneously 
forms  a  beam  in  the  desired-signal  direction  and  cancels  interference  in  a 
frequency-hopping  communication  system.  Each  weight-vector  iteration 
of  the  maximin  algorithm  is  specified  by  the  equations 

W(it+1)  =  W(A:)  +  -S-  *L1seW  _  gX#) 

A*)[  w  PnW  ' 


m= 


M 

W  ' 


_  ^  km 


p„(i) = i  £ 

m  i  =  {k-\)m  +  \ 


Ps(k)  = 


j  fc/n 

X  =  (fc5)m  +  1 


The  adaptation  constant  o;  usually  should  be  chosen  from  within  the 
range  defined  by 

0  <  a< 2  J  W(0)  |2pj  . 

In  these  equations. 


k  =  weight  iteration  number, 
i  =  sample  number, 

m  =  number  of  samples  between  weight  iterations, 

N  =  number  of  antennas  in  array, 

X(i)  =  input  vector  of  baseband  filter, 

3/r(0  =  real  part  of  output  of  baseband  filter, 
n(i)  =  input  vector  of  monitor  filter, 

V nR (0  =  real  part  of  output  of  monitor  filter,  and 

Pl  =  SINR  per  frequency  channel  at  an  antenna  input. 

If  the  fractional  bandwidths  exceeds  a  few  percent,  either  spectral  proc¬ 
essing  or  anticipative  processing  is  effective  as  a  means  of  frequency 
compensation.  Spectral  processing  divides  the  hopping  band  into  spectral 


regions  and  independently  adapts  the  weights  associated  with  each 
region.  The  anticipative  maximin  algorithm  for  the  auxiliary  weight 
vector  is 


w,('+i)=w«W- 


where  the  bracketed  factor  is  measured  by  the  auxiliary  adaptive  filter. 
The  maximin  algorithm  in  the  main  adaptive  filter  is  modified  by  requir¬ 
ing  that 


W(nk0)  =  W a(n  kj),  n  =  0, 1, ...  . 


In  these  equations, 

l  =  anticipative  weight  iteration  number, 
cq  =  anticipative  adaptation  constant, 
k0  =  number  of  weight  iterations  per  hop  in  the  main  filter, 
k-y  =  number  of  weight  iterations  per  hop  in  the  auxiliary  filter,  and 
n  =  hop  number. 
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